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ABSTRACT 

We show that smah sohds in low mass, turbulent protoplanetary disks collect into self-gravitating 
rings. Growth is faster than disk lifetimes and radial drift times for moderately strong turbulence, 
characterized by dimensionless diffusivities, ag < 10~^ — 10"'^ when particles are mm-sized. This range 
reflects a strong dependance on disk models. Growth is faster for higher particle surface densities. 
Lower gas densities and larger solids also give faster growth, as long as aerodynamic coupling is tight. 
In simple power law models, growth is slowest around 0.3 AU, where drag coupling is strongest 
for mm-sized solids. Growth is much faster close to the star where orbital times are short, with 
implications for in situ formation of short period extrasolar planets. Growth times also decrease toward 
the outer disk where lower gas densities allow greater particle settling. Beyond roughly Kuiper Belt 
distances however, solids are sufficiently decoupled from gas that dissipative gravitational instabilities 
are less effective. Turbulence not only slows growth, but also increases radial wavelengths. The 
initial solid mass in an unstable ring can be ~ lO^^M^ or greater, huge compared to kilometer sized 
planetesimals. Nonlinear fragmentation, which has not been studied in detail, will lower the final 
planetesimal mass. We consider applications to the asteroid belt and discuss the alternate hypothesis 
of collisional agglomeration. 



Subject headings: hydrodynamics — instabilities - 
and satellites: formation 



1. INTRODUCTION 

1.1. Summary of Paper I 

lYoudinl l|2005l hereafter Paper I) analyzed the linear, 
axisymmetric gravitational instability of solids (GIS, or 
GI for gravitational instabilities in general) subject to gas 
drag. Growth was described in terms of = ritgtop; the 
dimensionless measure of particle stopping time, tstop, 
where f2 is the orbital frequency; 
rameters: 



Qr = 



ttGY. 



and two stability pa- 

(1) 
(2) 



The Toomre parameter, Qt, measures the influence of 
particle velocity dispersions, c, in a Keplerian disk with 
particle surface density E, where G is the gravitational 
constant. The parameter Qn uses the particle sublayer 
thickness, /i, to measure the ratio of the Roche density, 
~ r2^/(7rG), to the particle space density, p = E/Zi. Pa- 
per I found that dissipation of angular momentum allows 
GIS for arbitrary values of the stability parameters, but 
growth is slower when they are large. 

Paper I also derived values of Qt and Qji generated 
by turbulent stirring of solids. Gas turbulence was char- 
acterized by a diffusive viscosity. 
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eddy turnover time, tQ. Unless otherwise noted, we take 
to — 1/0, i.e. orbital turnover times. 

The main equations from Paper I used in this work are 
briefly summarized. Our general dispersion relation (eq. 
[13] of Paper I, hereafter eq. [1.13]) defines the dimen- 
sionless growth rate 7 (eq. [1. 11]). Wavenumbers, fcf, 
(or wavelengths. At) of the fastest growing modes satisfy 
equation (1.17). Particle scale heights are determined by 
balancing vertical settling and turbulent diffusion (eq. 
[1.29]). The velocity dispersion (eq. [1.30]) includes di- 
rect kicks from turbulent fluctuations (eq. [1.25]]) and a 
contribution from epicyclic motion (eq. [1.32]). The re- 
sulting values of the stability parameters are plotted in 
Figure 1.3 (i.e. Fig. 3 of Paper I), and limiting analytic 
expressions for Qn and Qt are given in equation (1.36). 



1.2. This Work 

Since growth rates are always positive, our goal is 
to determine the conditions that give sufficiently fast 
growth times: 



4 



(4) 



(3) 



where Cg is the gas sound speed. The dimensionless dif- 
fusivity, ag, is defined by analogy with the accretion disk 
literature, but it refers to diffusion of material, not angu- 
lar momentum (see H5. 3. 31 for a detailed discussion). The 
analysis requires an assumption about the characteristic 



Growth should be faster than the 1-10 Myr lifetimes 
of gas rich T Tauri disks (Armitagc ct al. 2003), par- 
ticularly if solid cores are to accrete massive gas atmo- 
spheres. Since isotopic dati ng gives a 1 -3 Myr spread in 
age of meteoritic inclusions l|Amelin et al. 2002), there is 
no indication that planetesimals must form much faster 
than this. Second, growth shoul d be faster than aero - 
dynamic radial drift t imescales l|Weidenschillingl ri977al: 
iNakagawa et al.l 119861) . Also, the wavelength of unsta- 
ble modes must be smaller than the disk radius, R. We 
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TABLE 1 
Symbols 



— ; — ; — 

Symbol 


— — 

Reference Eq. 


— ^ 

Meaning 


fl t u 
'^^5 ''orb 


''orb ^ /I / J t 


orbital frequency, time 


/ 

''grow ? / 




growth time, dimensionless rate 


i-stop ) ' s 


GCJS. ^^^^^^^ii 


L/dl LiLlt; & LOU LJIllli Lllilt; 


''dr 




timescale for inward particle drift 


h, C 


eqs. (1.29,1.30) 


particle scale height, random speed 






stability parameters for solids 


ULg 


cq. 131 


turbulent diffusion parameter 


s 




particle surface density 


n 
H 


E//i 


particle space density 


a, ps 




particle size, internal density 


fp^ fg 




particle, gas enhancement factors 


h 


fcf = 27r/Af 


fastest growing wavenumber 


Af 


eq. 0^9) 


fastest growing wavelength 


Ag 


eq. |HJ 


traditional GI wavelength 




27rSAffl 


mass of solids in unstalale ring 




eq. ifel 


pressure parameter 






cylindrical disk radius, in AU 


Pg, Sg 




gas space, surface density 


Ms 




midplane gas fraction 


Amfp 




gas mean free path 


Qsonic 




"mixed" Toomre parameter 




-4-2 2 

log(Ts) 



implement these conditions as: 

igrow < O.lidisk = 10^ yrs , (5a) 

tgrow < tdr = R/ Vdr , (5b) 

Xf<R/2, (5c) 

where the timescale for inward drift of solids, tdr, is given 
in equation These conditions may not be the only 

ones that allow planetesimal formation, but they describe 
the limits of validity of our local analysis of GIS. 

This paper is organized as follows. We present gen- 
eral, model independent results in SJ2 We define simple 
powerlaw disk models in Growth for fixed Ug values 
is investigated in ^ The strongest levels of turbulence 
which allow GIS according to the above criteria (eqs. [SJ 
are explored in |S| Comparison to expected levels of tur- 
bulence is made in H5.8I The results are discussed in ^ 
Appendix 1X1 investigates the effect of fast eddy turnover 
times. 

2. MODEL INDEPENDENT RESULTS 

We first describe the instability in terms of ag and 
without assuming a specific disk model. We must specify 
a third parameter, Qsonic (eq. El)-^ We assume that 
eddies have orbital turnover times, but faster eddies are 
considered in Appendix^ This analysis provides a quick 
understanding of the general behavior of unstable modes, 
but does not allow evaluation of the criteria in equations 

2.1. Growth Times 

Figure ^ plots growth times relative to the local or- 
bital time, tgrow/^orb = 'iiT / "f . Thc black curves hold ag 
fixed. The stability parameters Qn and Qt for this case 
were plotted in Figure 1.3 {top). Growth is fastest for 

^ A two parameter description is possible with orbital turnover 
times since y/a^Qsonic always appears together. We opt for the 
familiarity of "a parameters" over a more compact description. 



Fig. 1. — Growth time for gravitational collapse of solids vs. 
stopping time. Black curves: Turbulent diffusivity, ctg, held fixed. 
Grey curves: Qr, inversely proportional to particle space density, 
fixed. Both curv es are for orbital turnover times and Qsonic = 
5 X 10^. See JTT\ 



weak turbulent diffusion (small ag) and marginal cou- 
pling, Tg « 1. Tight coupling retards collapse by slow- 
ing the terminal velocity of solids. Loose coupling low- 
ers the dissipation rate of angular momentum needed for 
collapse of long wavelength modes. We emphasize that 
growth times are always finite. 

Growth is also faster for higher particle surface den- 
sities, as self-gravity is stronger. Increasing S decreases 
growth times with a dependence between igrow cx: l/E^ 
and l/E for waves long and short (respecitively) com- 
pared to h. This follows from thc limiting analytic ex- 
pressions in Table 121 since Qsonic oc 1/E. These expres- 
sions, derived from the results of Paper I (using eqs. [1.18, 
1.19, 1.21, 1.B2a]), explain the general features of the con- 
stant ag curves. The one exception is the rapid dynam- 
ical collapse for ag = 10~^ and Ts > 1. In this case 
Qb. ^ Qt ^ 1 and traditional GI dominates dissipative 
growth (see eq. [I-Bl] and surrounding discussion). 

Grey curves (in Fig. ^ plot growth at fixed Qu. 
Growth in under 10*^ orbits is possible for Qji < 800, 
densities nearly a thousand times smaller than the Roche 
density. Compared to the standard Qr < 1 criteria, dis- 
sipation allows GIS for much stronger turbulence. 

2.2. Radial Wavelengths 

Figure El (bottom) compares the most unstable wave- 
length, Af — 27r/fcf to the standard (non-dissipative) 
wavelength of GI, 

Ag = 2TT^GY./n^. (6) 

We find Af S> Ag unless turbulence is very weak or (con- 
sulting Tab. El S is very large. Such long wavelengths 
are allowed by angular momentum dissipation. The ini- 
tial ring mass ~ 27rI]Afi? vastly exceeds the traditional 
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TABLE 2 

Limiting properties of unstable modes 



Quantity 




Regime'' 




Ts < 1, kfh < 1 


Ts < 1, kfh > 1 


> 1, fcfh < 1 


^igrow 






°9Qsonic^'' 


Af/Ac 




2/3„4/3 -1/3 


"9Qsonic 


kfh 


Kr,)-V2Q-1.^ 




Kr.)-l/2Q-l.^ 



^For orbital turnover times, tq = 1. 

''The Ts 3> 1, fcf/j <SC 1 case is not particularly relevant and was omitted. 




Ji 102 



Fig. 2. — Top: fastest growing wavenumbcr times scale height 
vs. stopping time for fixed ag. Bottom: most unstable wavelength 
compared to the standard GI wavelength. Both plots assume or- 



Figure|21 {top) plots kfh, measuring the size of waves 
relative to the layer thickness.^ Since kih < 2tt, our 
vertically integrated model is a good approximation. Our 
limiting analytic formulae (Tab. ^ apply for kfh <^ 1 or 
3> 1. Figure 13 indicates which case is more applicable. 
When kfh ~ 1, behavior is intermediate and requires 
numerical evaluation. 



in 




0.1 



1.0 10.0 
R [AU] 



100.0 



Fig. 3. — Dimensionless stopping time vs. disk radius for mm 
(black) and cm (grey) sized solids with a material density of 3 
gm/cm^ in a standard gas disk (solid lines) and with gas depleted 
(dashed lines). Marginal coupling is indicated (dotted line). Kinks 
demark the transition between Stokes (small R) and Epstein (larger 
R) drag, which in reality would be smoothed. 



estimate ^ SAf cor respon ding to km-sized planetesimals 
ijGoldreich fc Wa rd ll97.lD , but subsequent fragmentation 
is likely (see M5.2|) . The weak dependence on bodes well 
for studies that will include a dispersion of particle sizes. 
A detailed analysis is needed, but a range of particle sizes 
could collect in the same annulus with larger solids (i.e. 
closer to Tg = 1) collapsing faster. 



3. PROTOPLANETARY DISK MODEL 

Our numer ical analysis uses a standard mini mum mass 
model ( Wcidcnschilling 1977b; Havashi 1981; our nota- 
tion follows Youdin fc Shu.2002. hereafter YS02) for the 

^ Since Af/Ag = Q n/ (kfh), the two plots in Figurcl2l differ only 
by a factor of Qr. 
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particle and gas surface densities and the sound speed: 



S = /pl0tn-3/2 g/cm^ 
E<, = /3l700n7-3/2 g/cni3 
Cg — XvD^^I'^ km/s 



(7) 
(8) 
(9) 



where vo = R/PM gives stellocentric distance in AU, 
and fp and fg are enhancement (or depletion) factors 
for solids and gas, respectively. Our reference model has 
fg = fp = l. We treat fp and fg as constants to under- 
stand the behavior of simple power law models. How- 
ever many enrichment or depletion mechanisms deviate 
signficantly from power law behavior, including: ice con- 
densation (many models give fp — ^ outside the snow- 
line), particle pile-ups generated by inw ard radial drift of 
solids rYS02: Wondin fc Chiang"2004'. hereafter YC04), 
and photoevaporation of gas (ThrooD & Bally 200^. 

Particles are assigned a uniform size a (1 mm unless 
stated otherwise) and density Ps — i g/cm^. These are 
typical properties of the abundant meteoritic inclusions 
called chondrules. It is a realistic concern that a disper- 
sion of particle sizes might yield less efficient collapse. 
We plan to study this issue in the future , but n ote that 
chondrules are strongly size-sor ted (TavloJ200l p. 117). 

For a < (9/4)Amfp, where 



A 



mfp 



IfZ^uj^-''' cm. 



(10) 



is the gas mean free path, the stopping time is given by 
Epstein's law: 

rfP = V2^p,a/Eg « 4 x l{)-\,,,^w'"^ / fg, (11) 

where amm = a/(l mm). Unless stated otherwise, 
numerical estimates will use Epstein's law. For a > 
(9/4)Amfp, i.e. small Knudsen numbers. Stokes' drag law 
applies and 



St 



'2tt- 



9E„A 



9 X 10 



g -^mfp 



mm 



V0.3/ 



-5/4 



(12) 



independent of gas density. Figure |21 plots the stopping 
time for both mm and cm-sizes. The effect of gas deple- 
tion (to fg = 0.1) is demonstrated. Since rf* decreases 
with disk radius, while r^P increases, a given particle is 
most tightly coupled at the transition between the two 
regimes. We do not consider particles large enough for 
turbulent drag laws to be relevant. 

Other relevant quantities take the following values 

Q,o„ic = Cgf]/(7rGS) « 9.5 x lO^tu-i/V/p , (13) 



Qg=Cg^/{^GY.g) « bQVJ-^'^/fg, 



dP/dR 

^ 2n^Rpg 

1 + ^9^' . 



1.8 X 10 



-3 



yrs, 



(14) 
(15) 

(16) 



where pg — Pg/{pg + p) is the midplane gas fraction. 
Large Qg is consistent with our assumption that the gas 
disk as a whole is not gravitationally unstable. The nu- 
merical expression for idr assumes tight coupling and Ep- 
stein drag. Drift times are longer than nominal values if 
particle inertia dominates in the midplane to give fig < 1. 
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Fig. 4. — Growth times vs. 
erence model [fg = fp = a^n 



radius for fixed an 



1) for Og 



10 



Top: our ref- 
(solid lines), 



10"" (dashed), and 10"' (dot-dashed). Sharp maxima denote a 
transition between drag laws. Bottom: with Qg = 10~^, the disk 
model is varied by: adding solids (dashed line), depleting gas (dot- 
ted line), increasing particle size to 1 cm (dot-dashed line), and all 
of the above changes (triple-dot-dashed line). The fg = 0.1 model 
overlaps with the a = I cm curve for ro > 0.8 and the reference 
model for zu < 0.2, because Ts is identical in these regimes. 



4. GROWTH AT FIXED Og 

Figure 21 (top) plots the growth times in our reference 
model for several values of ag. The general results that 
growth is faster for weak turbulence and for near unity 
are confirmed. The sharp peak at Ci7 = R/AU « 0.3 
occurs at the transition between Epstein's and Stokes' 
laws, where is minimized. 

The detailed radial dependence is explained by the for- 
mulae in Table 121 for <C 1: 



OigQ 



2 

sonic 



-I7T73/2 



/(27r) 



-3/2 3/2 



7(4^) 



if kih < 1 
if kih > 1 



(17) 

Close to the star, in the Stokes regime, growth times de- 
crease steeply, roughly as tgrow oc R^ . This is a com- 
bined effect of faster orbital times and shorter Amfp, 
which weakens viscous drag coupling. Moving radially 
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Fig. 5. — Growth times vs. drift times for the same models as 
Figure El 



outward from the drag transition, where Epstein drag 
applies, growth times also decrease. Slower orbital times 
are compensated by the decrease in Ts at lower gas densi- 
ties. The radial decrease in Qsonic, a result of our choice 
of power laws, also favors growth at large radii. When 
Ts > 1 at sufficiently large radii, zu > 180(/g/amm)^^'^, 
looser coupling gives longer growth times. Even if the 
disk is not truncated, rapid GIS will not occur at very 
large stellocentric distances. 

In Figure 01 (6ottom), as = 10~^ is fixed to show the 
effects of other model parameters on igrow Enhancing 
the particle surface density via fp (dashed curves) gives 
stronger self-gravity and faster growth. Increasing par- 
ticle sizes, a (dot-dashed curves), or decreasing the gas 
surface density via fg (dotted curves) gives larger Tg.^ If 
T, < 1 these changes shorten growth times, but once 
Ts > 1, the same changes give slower growth. This 
explains the increase in tgmw beyond ~ 30 AU in the 
fg — 0.1 and a — I cm curves (and beyond ~ 8 AU when 
both changes are made). When multiple effects favor- 
able to GIS are included (gas depletion, larger particles, 
and/or particle enrichment) growth times are quite short 

^ Except T^* is independent of gas density. 
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fQ = 0.1 

f ^"^^ 
reference 



0.1 



1.0 10.0 100 
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Fig. 6. — Maximum levels of turbulent diffusion for which GI 
is relevant vs. disk radius for several disk models differentiated by 
linestyle. Line colors indicate the limiting factor, red for radial 
drift, blue for disk li fetim es, and green for wavelengths comparable 
to disk radius. See H5.ll 



in the inner disk, as seen in the curve labeled "all." 

Figure compares growth and drift times for the same 
models as Figure ^ In many cases igmw > ^dr, a viola- 
tion of condition l|5b(l , though many of these modes have 
igrow > 10^ yrs as well. The next section shows that drift 
rates do often set the most stringent constraint on ag. 
For now, note that the effects which g ive faster igj-ow also 
decrease tgrow/^dr- This is not immediately obvious since 
growth and drift times both decrease toward marginal 
coupling, but tgi-ow does decrease faster. Thus conditions 
which give faster growth also help growth outpace drift. 

5. ALLOWED LEVELS OF TURBULENCE 

5.1. Calculation o/amax 

We now calculate amax, the largest value of Og that 
satisfies all of conditions in equations Figure |B| plots 
Q!max against disk radius for the same models considered 
in Figure 01 plus a high enhancement, fp ~ 20, case.'* 



At 1 AU, 10"'^ < Qfmax < 10"^, showing that allowed 
levels of turbulence can be strong, but are model depen- 
dent. Drift is often the limiting factor that determines 
Q^max, as indicated by the line colors. Since particle pile- 
ups naturally enhance the surface density of solids on tdr 
(YS02, YC04), they could act in concert with GIS. A 
non-local model is required to explore this possibility in 
detail. Drift also introduces shorter wavelength instabil- 
ities which might co ncentr ate particles faster than t^ r 
ijGoodman PindoH l20?i?it lYoudin Goodma,nl 1200^. 
but their effect is not seen in the current single fluid, 
2D analysis. 

* The fp = 20 case is not well motivated in the outer disk where 
particle pile-ups have less of an influence. Also the large total mass 
of solids would require efficient ejection from our solar system. 
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Fig. 7. — Properties of unstable modes vs. radius for the strongest 
allowed turbulence. Linestyles and colors have the same meaning 
as Figure|S] Growth times (top), ratio of wavelength to disk radius 
(middle), and mass of unstable modes in Earth masses (bottom) are 
all maximum values that would be smaller for weaker turbulence. 
The reference, fp = 4, and fp = 20 models overlap on the top plot. 



The disk lifetime (condition I5a|l can set ag in two 
cases. The first is near the Epstein-Stokes transition 
where stopping times are particularly short. The second 
is in the outer regions (i? > 40 AU) of gas depleted disks, 
where > 1. The restriction that waves are shorter than 
the disk (condition |2c|l is relevant if conditions are very 
favorable to GIS, and amax is large. For instance, the 
model "all" has wavelengths comparable to R between 
- 2—10 AU. 

5.2. Properties of "Slowest Growing Modes" 

Figure H plots properties of unstable modes vs. R for 
ag = ttmax- This gives upper limits on growth times, 
wavelengths, and masses of unstable annuli, all of which 
would be smaller for weaker turbulence. The top plot 
shows that growth times can be as long as ~ 10^ yrs for 
a = 1 mm, but must be < 10^ yrs for a — 1 cm.^ 

^ Note that t^,. (which equals tgrow for the red [black in print 
version] curves) is not decreasing ex as one expects for an 

individual particle. This is due to greater particle settling and 



Collision times could be shorter, but this does not, by 
itself, qualify as a weakness of the GI hypothesis when 
compared to collisional agglomeration. First, we empha- 
size that we are quoting the longest allowed growth times. 
Second, the collision time for large bodies is quite slow, 
especially in the outer disk: 



^coll — 



ApsO 



10"— n7^ yrs, 
km 



(18) 



where the numerical estimate uses c « fJ/i, appropri- 
ate for loose coupling (§1.3.3). Collision times will be 
shorter if a populatio n of small bodies is maintained 
ijGoldreich et alJl2004j) . but this requires fragmentation. 
This brings us to the important point that collisional 
growth times involve a balance between fragmentation 
and agglomeration. With fragmentation included, it is 
not clear that the collisional growth time to km sizes is 
finite, let alone faster than fjr (very rapid in the 10 cm — 
1 m range) or disk lifetimes. The collisional hypothesis 
is discussed further in fJB| 

The middle plot shows that wavelengths are not negli- 
gible compared to R. This raises the possibility of high 
resolution mm-wave observations, perhaps by ALMA. 
Prospects would be particularly good for gas depleted 
or particle enhanced disks around ^ 30 AU. For non- 
dissipative GI of solids, \g/R ~ 10~^/pCC7^/^, and there 
was little hope of detection. 

Upper limits on the mass of unstable modes are plot- 
ted in Figure (ftottom). These masses include a com- 
plete annulus, M\ = 2ttRX['S,, because with t. 



grow 



> 



iorb initial growth must be axisymmetric. The reader 
should not be too concerned that these masses dwarf 
the "canonical" values for kilometer-size planetesimals, 
M\ ~ SA^, obtained from no n-dissipative stability cri- 
teria iColdreich fc WardllTOTl . 

Non-linear effects should give subsequent fragmenta- 
tion of the initial annulus. As noted by Ward (1976), 
when "the local surface density increases, material be- 
comes more susceptible to shorter-scale local dynami- 
cal instabilities." Indeed the ring should fragment az- 
imuthally, and perhaps radially, when tgrow ^ ^orb • Frag- 
mentation will be enhanced by mass loading, the feed- 
back of particle inertia on the gas turbulence (see tl5.3.3|l , 
and the increasing rate of inelastic collisions. As in star 
formation theory, the loss of support should lead to hier- 
archical fragmentation. Without a detailed understand- 
ing of these processes, we cannot estimate the planetesi- 
mal masses produced by GIS. 

Nevertheless, by locking large amounts of solids in a 
mode for many orbital times, GIS would help preserve 
the chemical gradients in the solar nebula. If drift com- 
presses relative separations, as one would expect from 
particle pile-ups, then gradients could be amplified as 
well. Implications for the asteroid belt are discussed in 

m 

5.3. Can Turbulence Prevent Gravitational Instability? 

We now compare our values of Q;,nax, which range from 
10^^ to 10~^ depending on disk model and radial posi- 
tion, to the expected levels of turbulence in protoplane- 
tary disk midplanes. We first recall that previous analy- 

smaller fj.g at large R. 
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ses (e.g. lWeidenscliillindlT995j) using the Qr < 1 stabil- 
ity criteria placed stronger constraints, 

«9 < 5 X 10-i2tn2a„„j2//s , (19) 

for tight coupling and Epstein's law. Thus GIS is relevant 
for turbulence stronger by at least 5 orders of magnitude 
when dissipation is included. 

5.3.1. Turbulence Driven by Vertical Shear 

GIS involves particles settling to the disk midplane, 
and turbulence arising from this stratification is well 
studied. The source of instability is vertical shear in 
the azimuthal velocity. Particle-free gas has an orbital 
speed that is slower than Keplerian by "qvK due to ra- 
dial pressure support. As the gas is loaded with (tightly 
coupled) solids its orbital speed increases toward the Ke- 
plerian value. Thus a vertical density gradient drives a 
velocity shear which can trigger turbulence as with the 
Kelvin-Helmholtz instability. Semi-analytic analyses of 
these shear layers d erive a thick ness h k, r]R/2 in the 
tight couphng limit l|Seki vail 19981 YS02). For looser cou- 
pling, one expects stronger tu rbulence and thinn er par- 
ticle layers, in agreement with ICuzzi 

YS02 showed that only a critical surface density, 
Ecrit ~ 2pg/i, of small solids can be stirred. Any ex- 
cess solids will form an overdense midplane layer, which 
would be subject to rapid GIS. Since the current work 
(and Paper I) does not include mass loading effects, the 
saturation effect does not occur. Thus our growth rates 
are probably too conservative, particularly for solid to 
gas ratios enhanced relative to cosmic abundances. 

In light of the current work, we ask whether stirred 
solids might also participate in GIS. With h k, riR/2, 
the parameter Qji « 250/ f p. The constant Qr curves of 
Figure^show that even with /p = 1 marginally interest- 
ing growth is possible if Tg > 10~^. With some particle 
enhancement, and a lower Qn, growth would be faster 
and relevant for more tightly coupled solids. Thus GIS 
may take hold even if E < Ecrit- Ongoi ng effor ts to un- 
derstand vertical shear i nstabilities ll|Garaud fc Lin 2004; 
IGomez fc Ostrikeil2005l) are improving understanding of 
this stirring mechanism and its effects on GIS. 

5.3.2. Accretion Related Turbulence 

More generally, protostellar disks are thought to be 
turbulent because mass accretion requires an angular mo- 
mentum transport mechanism that is thought to be tur- 
bulent in nature. The observed accretion rate of classical 
T-Tauri stars (CTTS) at 1 Myr is M - IQ-^ Mo/yr with 

order of magnitude scatter and a trend for M to decrease 
with age ijHartmann et alJ ll998^. These accretion rates 
are reproduced by vi scous accre tion disk models with 
10-3 < aL < 10-1 ()D'Alessio et all 11999ft . where the 
diffusivity of angular momentum, aL, may differ from 
ag (see M5. 3.3(1 . The relation between M (diagnosed ob- 
servationally with UV veiling and emission line profiles) 
and aL (from theoretical modelling) is not unique. Thus 
the decrease in M with age could be due to weaker tur- 
bulence, lower disk masses, or a combination. 

The magneto-rotational instability is a leading candi- 
date for the driver of angular momentum transport mech- 
anism and disk accretion. (Torques from density waves 
in massive gas disks are another possibility.) Ideal MHD 



simulations indicate that the MRI generates 5 x IQ-"^ < 
aL < 0.5 ijStone et al.ll200d() . Much lower values are 
possible, including complete quenching, when resistive 
effects are important, i.e. the ionization fraction is too 
low. In the "layered accretion" scenario accretion flows 
are confined to ionized surface layers, leaving t he resistive 
midplane relatively quiescent l|Gammielll996j) . 

In summary there is observational and theoretical sup- 
port for characteristic values of 10-^ ^ ocl ^ 10~^ for 
angular momentum diffusion, but lower values are al- 
lowed. We now consider the relation between aL and 
diffusion coefficients relevant to GIS. 

5.3.3. Relation Between Diffusivities 

For GIS we are mainly interested in the diffusivity of 
solids, as- Since particles respond to turbulent fluctua- 
tions of gas, we expressed our results in terms of ag, using 
a Schmidt number Sc = ag/as- The relation between ag 
and aL, needed for comparison to stellar accretion rates, 
is characterized by a Prandtl number, Pr = aL/ctg. Ad- 
ditional relations are possible for anisotropic diffusion, to 
c ompare radial and vertic al transport. 

IGarballido et al.l l|2005|) compared the radial diffusion 
of angular momentum and a passive contaminant (gas) 
in MHD turbulence and found Pr ~ 10. This result 
is readily explained by Maxwell stresses, which trans- 
port angular momentum without generating large ve- 
locity fluctuations. This finding favors GIS. However 
IJohans en & Klahr (2005) included vertical transport and 
particle coupling to measure a^/as = Pr Sc < 1. These 
results appear to be a contradiction if one expects Sc > 1, 
i.e. particles to be less diffus i ve tha n gas. 

We used the |C UZZl et all (|1993|) result for Sc, equa- 
tion (1.26), which is greater than unity due to the imper- 
fect coupling of particles to eddies. If particles prefer- 
entia lly concentrate in quiescent nodes of the turbulent 
flow (|C uzzi et aT]l2001j) , then diffusion of solids would be 
less efficient, and Sc would be even larger. 

Mass loading, the feedback of particle inertia on tur- 
bulence, should also decrease stirring efficiency. Mass 
loading should become significant for p > Pg (YS02), or 
Qr < 140ro-i/4//g. While GIS is already significant 
at these densities without mass loading (see Fig. QJ, it 
is likely more rapid. Mass loading will certainly become 
significant as GIS develops, and the loss of turbulent sup- 
port could lead to hierarchical fragmentation (see H5.2|l . 

To conclude, a variety of effects can decrease the effi- 
ciency of particle stirring in turbulent disks. These in- 
clude mass loading. Maxwell stresses, layered accretion, 
and turbulent concentration. The role of these effects 
is only beginning to become clear. Even with the pes- 
simistic assumption that all the above effects are irrel- 
evant and aL = ag, disks must only be somewhat less 
active than characteristic CTTS values for GIS to be rel- 
evant. 

6. DISCUSSION 

We investigated the gravitational instability of solids 
(GIS) in protoplanetary disks with dissipation provided 
by gas drag including the stirring of particles by turbu- 
lence. GIS is faster for weak turbulence, high abundances 
of solids, and particle stopping times comparable to or- 
bital times. We conclude that GIS is a viable mechanism 
to form planetesimals in disks that are not too strongly 
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turbulent. More importantly, as knowledge of physical 
conditions in disks improves through observation, the- 
ory, and simulation; the framework developed here can 
make firmer predictions about the role of GIS. Possible 
refinements to our dynamical model were discussed at 
the end of Paper I. 

An alternate hypothesis is pl anetesimals form by col- 
lision al agglomeration alone ijWeidenschilling fc Cuzzil 

I1993D. The physic s of dust coagulation 

ijDominik fc TielensI Il997j) is sufficiently complex 
that predictions for growth beyond even fun sizes are 
difficult. Existing theories do not include a realistic 
balance between constructive and erosive impacts over 
the relevant size range, from /im — km. Admittedly this 
a daunting task, and agglomeration cannot be ruled out 
simply because it computat ionally i nt racta. ble. However, 
as argued in YS02 and lYoudiiJ l)2004|) . laboratory 
experiments and basic physical considerations do not 
give much reason for optimism. Various mechanisms 
have been proposed to promote coUisional growth. 
For instance, gas flow through a porous body might 
return coU isional fragments fo r repeat impacts a t lowe r 
velocities l|Wurm et al.l 12001: iSekiva fc Takedal l2n(il . 
Without a broader theory of coagulation in protoplan- 
etary disks, the relevance of such sticking enhancement 
mechanisms is not clear. 

Observational evidence for grain growth is stronger. 
The most robust explanation for sub-mm spectral in- 
dices of protostellar d isks is particle gr owth to sizes 
> 3 mm (I Calvet et al.ll2002HChi ang"2004: 'Draine'l 20n5D . 
This growth could be from coagulation of dry and/or icy 
grains, but might also involve partial melting, as in the 
format ion of refractory meteoritic inclusions li ke chon- 
drules ^Shu et a,l .1120011 IDesch fc CoTmo]lvll200^. What- 
ever the mechanism, growth to mm sizes aids settling 
to the midplane and abets GIS. 

We found that GIS develops as rings with radial wave- 
lengths much longer than for traditional GI. The large 
initial solid masses of unstable annuli does not neces- 
sarily imply large planetesimals would be formed. Be- 
cause subsequent fragmentation is not yet understood, 
we do not predict final planetesimal masses. The 
potential combination of long growth times and AU- 
scale wavelengths increases the prospects for observation. 
ALMA might image giant planets in the gaps they cre- 
ate (|Wolf fc D'Aneclo 2005). Detecting a ring of solids 
would be more difficult, since sharp edges are not likely 
with a distribution of particle sizes. 

The implications of GIS for the asteroid belt should 
be considered. Known meteorite falls from the asteroid 
belt (numbering > 22, 000) are associated with a rela- 
tively small number, 100 — 150, of distinct "parent bod- 
ies." Some of these correspond to well-known asteroid 
classes. It has been noted that the term "parent body" 



could be misleading if many asteroids form from nearly 
identical material fsee iBurbine et all 120021 for a review 
of meteoritic parent bodies). The chemical homogene- 
ity within a given group of asteroids might be due to 
> lO** years of mixing in a slowly collapsing ring prior to 
non-linear fragmentation into many planetesimals. 

The strong chemical zo nation between classes in the 
asteroid belt l)Tavloill2001l p. 280) might reflect the com- 
bined action of GIS, which segregates distinct annuli, and 
drift pile-ups (YS02, YC04). By bringing neighboring 
rings closer together, pile-ups could steepen pre-existing 
compositional gradients. These possibilities are rather 
speculative, but serve as a reminder that, despite 4.5 Gyr 
of coUisional and dynamical evolution, the arrested de- 
velopment captured in asteroid belt provides invaluable 
clues for and tests of planetesimal formation theories. 

Further circumstantial support connecting radial drift 
(but not necessarily GIS) to planetesimal formation is 
found in Jupiter's atmospheric abundances. The uniform 
enrichment of noble gases and nitrogen points to a low 
temperature, ~ 30 K, origin for Jupiter's planetesimals 
ijOwen et al.lll999() . 

The growth rates of GIS vary significantly with disk 
radius. Slowest growth occurs near Mercury's orbit for 
mm-sized solids in a standard disk model. The sharp 
decrease in growth times at closer distances to the star 
is encouraging for the in situ formation of hot Neptunes 
and Jupiters, esp ecially when th ere is a nearby stellar bi- 
nary companion l|Konacldl200,'iD . Growth times were also 
found to decrease toward the outer disk until decreasing 
gas density causes > 1, and growth times increase. 

This generally suggests that planetesimal formation 
at > 100 AU (for mm-sized solids) might be difficult. 
However it does not obviously predict that planetesimals 
should only form inside ~ 50 AU, as indicated by Kuiper 
Belt observations (Allen et al.,2002). By contrast, YS02 
and YC04 argued for a hard edge based on the pessimistic 
belief in threshold criteria for GI. While it's a welcome 
change to have a theory that threatens to overproduce 
planetesimals, several scenarios could produce an edge 
inside ^ 50 AU. Particles may migrate to the inner disk 
while turbulence is too strong for GIS. Alternatively an 
external O star may rapidly photoevaporate the outer 
disk. 
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APPENDIX 
A. FAST EDDIES 

In §1.3.4.2, we considered the stirring of particles by fast eddies (FE) with short turnover times. Figure EH plots 
growth times for this case, whose stability parameters were plotted in Figure 1.3 (bottom). We compare to the growth 
rates found for orbital turnover times (OTT) in Figure ^ For growth with Ug fixed (black curves), FE gives faster 
growth than OTT when > 1. This is easy to understand as both Qji and Qt are smaller for FE since loosely 
coupled particles respond poorly to rapid fluctuations. 

For Ts < 1, the growth is slower for FE than OTT. Here coupling to the faster eddies is strong enough that velocity 
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Fig. A8. — Similar to Figure ITl except for faster eddies. Linestyles have the same meaning as in Figure^ Growth is slower for tight 
coupling but faster for loose coupling. The drop in growth rates for Ts > 10 in the Oig = 10~® curve occurs because Qt drops below unity 
(see Fig. 1.3, bottom). 



dispersion, and Qt, are larger. The higher densities (lower Qfj) for FE is less relevant because Q'^/Qr ^ 1. This 
corresponds to the thin disk, kfh <C 1, regime where growth is controlled by the velocity dispersion, not the density 
(see §1.2.4). 

However these growth time estimates for FE are probably too conservative. We showed in §I.C that collisional 
damping can be ignored when c/ {Q.h) = Qt/Qr ^ 1, which is the case for OTT (when < 1). But when c/ (Vih) ^ 1, 
as is the case for FE (see §1.3.4.2), we showed higher collision rates give significant dissipation. This would give faster 
growth for FE. Thus a more detailed analysis is required to determine if FE impedes or promotes growth (relative to 
OTT with constant ag) for ^ 1, but it clearly promotes growth for Ts > 1. 

The growth times at fixed Qr (grey curves) are much longer for FE, simply because it takes much more vigorous 
turbulence (larger Ug and thus Qt) to stir solids to a given height with FE. This demonstrates again that dissipative 
GI cannot be understood with a single stability parameter. Note that the fluid approximation is questionable when 
igrow < istop, i-e. igrow/^orb < Tg / {2?:) , and a kinetic theory treatment would be mor e ap propriate. This technical 
comment applies to dynamical collapse in the lower right regions of both Figures HI and IXsl 
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